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1 Introduction 

The spacetime as we know it is four jiimensional. A fifth dimension was pos- 
tulated by Kaluza and Klein (KK) EJ in the 1920's, in an attempt to unify 
electromagnetism and gravity. If a fifth dimension exists, it must be very 
small or else it would have been seen. To estimate how small that is, recall 
that the allowed wave numbers in the fifth dimension are integral multiples of 
1/-R, the inverse radius of that extra dimension. Excited wave numbers are 
seen as massive particles in our four-dimensional world, with a mass gap 1/ R. 
Since no such KK particles are known up to about 1 TeV, the radius R must 
be smaller than 1 (TeV)^^ = 2 x 10^^^ m. This is the picture of the fifth 
dimension until a few years ago. 

Recently a different view emerged 13 which allows the extra dimensions to 
be large. Partly motivated by the discovery of D-branes in string theory 0, it 
postulates all Standard Model (SM) particles to be permanently confined to 
our four-dimensional world, the so-called '3-brane'. Being confined they can 
have no KK excitations, so the previous bound on R is no longer valid. On the 
other hand, SM singlets such as gravity and right-handed neutrinos are allowed 
to wander into the extra-dimensional world, 'the bulk'. As a consequence, the 
inverse-square law of gravitational force is modified at separations less than R. 
Such a deviation has not been detected down to submillimeter separations □, 
from which we obtain an upper bound of R to be about 10^^ m, not the much 
smaller 10~^^ m discussed before. Moreover, if there are n extra dimensions, 
the fundamental energy scale will change from Mp ^ 10^^ GeV to Mf = 
[M|,/(27ri?)"]i/("+2). For n = 2 and R = 0.4 mm, this can be as low as 1 
TeV, leading to many observable consequences. However, if n = 1, Mf = 
2.3 X 10* GeV remains quite beyond our present capability to reach, then 
neutrino physics is probably the only way to detect the extra dimension. This 
is why neutrino physics is so important and so interesting in this connection. 
In this talk I will discuss some of the consequences for neutrino oscillation when 
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there is only one large extra dimension present. The small extra dimensions 
decouple and will not be taken into account. For simplicity I will assume the 
fifth dimension to be flat, though a curved scenario is also interesting^. 

2 Neutrinos Are Different 

Neutrino is unique as an extra-dimensional probe because the right-handed 
neutrino is a SM singlet. A large radius R in the extra dimension creates a 
small energy scale 1 / R that only neutrinos can see, which may be why neutrinos 
have such exceptionally small masses. For R = 0.4 mm, the energy 1/R is half 
a millivolt. Whether this is really the natural energy scale is unimportant for 
what follows because I will treat the neutrino masses as parameters. 

A right-handed neutrino roaming in the bulk is derived from a 5-dimensional 
Dirac field, and is sterile. It gives rise to a KK tower of left-handed and right- 
handed neutrinos. The minimum content of a five-dimensional theory consists 
of three active brane neutrinos and a sterile bulk neutrino. Depending on the 
model, there may of course be more brane neutrinos which are sterile, and/or 
more bulk neutrinos. 

Quark mixing via the CKM matrix is small, but at least some of the mix- 
ings for active neutrinos is large. One might be able to expalin this difference 
if the additional amount of neutrino mixing is derived from their coupling to 
the bulk neutrino(s). 

In the presence of a bulk neutrino in a large extra dimension, neutrino 
oscillations possess special characteristics. The infinite number of states in the 
KK tower generally leads to a very complicated oscillation pattern. Besides, 
these neutrinos are sterile so they have neither charge nor neutral current in- 
teractions with matter. These are features detectable from neutrino oscillation 
experiments. 

I will discuss two simple models IB in some detail to illustrate these char- 
acteristics. These models are instructive though possibly too simple to be 
realistic. I will also discuss more complicated models later in the section on 
phenomenology. 

3 Mixings and Oscillations 

Before discussing these models let me first review the usual formalism of 
neutrino oscillations. Let Vejiy^^iyr be the active neutrinos in the flavor ba- 
sis, namely, those obtained directly from the charged current decays of W^. 
Through their mutual interaction and perhaps interactions with sterile neutri- 
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nos that may be present, they mix to form the mass eigenstates vy. 



A 



(1) 



There are as many eigenvalues A as there are total number of neutrinos, active 
and sterile ones both counted. The sum is taken over all these eigenvalues, 
which consists of 3 + s terms in a four-dimensional world with 3 active and 
s sterile neutrinos. In a five-dimensional world, there are always an infinite 
number of sterile neutrinos so it is an infinite sum. 

The transition probability into brane species «, after the incoming neutrino 
of species j and energy E has traversed a distance L, is equal to Vij (r) = 
|^y(T)p, where t = L/2EE?. I will use \/R as the basic energy unit and 
keep all the other parameters dimensionless. The transition amplitude is given 
by 



(2) 



so the transition probability is 



P.,(r) = l-2 



A">A' 



2 Re( U*yU,yU^y'U*y, 1 sin^ f ^AA^ r 



lm(u*yUjX'U^y,U*y, ) sin(AA2 r) 



(3) 



where AA^ = A"^ — A'^. It is clear from these formulas that the oscillation pat- 
tern gets more and more complicated as the number of eigenvalues A increases. 
In the presence of an extra dimension, this number is always infinite, so the 
oscillation pattern is very complcated indeed. An exception occurs when the 
widths of the KK resonances are large compared to their separations, in which 
case the resonances all merge into a continuum background and the oscillation 
pattern becomes relatively simple again. This situation will be discussed in 
Sec. 4.4. 



4 Two Simple Models 

I will discuss two simple models, DDci and ADDM0, slightly generalized from 
their original forms. Both are assumed to contain one Dirac bulk neutrino, as 
well as / left-handed brane neutrinos coupled to the bulk neutrino via Dirac 
mass terms proportional to di/R. These two models differ in that lepton 
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number is conserved in the ADDM model, but is violated by the presence of 
Majorana masses rrii/ R in the DDG model. The parameters d^, rrii are assumed 
to be real but otherwise completelv arbitrary. In units of 1/i?, the neutrino 
mass matrix in the ADDM model Q is 
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(4) 



The rows and columns of this matrix are labelled respectively by the left- 
handed and the right-handed neutrinos. The brane neutrinos are purely left- 
handed, but the KK tower of bulk neutrino contains both left-handed and 
right-handed components. They have masses n/R, with n = 0, 1, 2, • • •. The 
n = mode is special, because its left-handed component is decoupled from 
everything else, so it does not appear in the mass matrix. Its right-handed 
component does couple to the / brane neutrinos and it occupies the first column 
of the matrix. In short, the columns are labelled by the n > modes but the 
rows are labelled by the / brane neutrinos- followed by the n > modes. 

The brane neutrinos of the DDG modelQ are Majorana so rows and columns 
of its mass matrix are both labelled by all the brane neutrinos and all the modes 
of the bulk neutrinos. In imits of l/R, the mass matrix of this model is 
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Note that the nth KK mode of the ADDM model is a linear combination of 
the ±n modes of the DDG model. Only one combination occurs in the mass 
matrix M' because the other combination decouples. 
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We shall assume the 2/ parameters rrii and di to be real, in which case 
the mass matrix M in the DDG model is real and symmetrical. The unitary 
mixing matrix U (actually real orthogonal in this case), and the eigenmasses 
A/i?, can both be obtained by diagonalizing M. The resulting eigenvalues A 
satisfy the characteristic equation 

itan(^A) = f]-^, (6) 

i—1 

whose graphical solution is illustrated in Fig. 1 for / — 3 and two different 
values of d^ = X]f=i '^l- We shall denote tan(7rA)/7r by 77(A), and the sum on 
the right-hand side of (||) by d'^r{X). We shall also write di = dei. Therefore 
E,e2 = landr(A)=E.e|/(A-™,).^ 

The components of the corresponding eigenvector 

i= {wi,W2,---,Wf,Vo,Vi,V-i,V2,V-2r--)'^ (7) 

are 

m,(A) =d,/(A-m,) (i = 1, 2, • • • , /), 

^;„(A) = ry(A)/(A - n) (n e Z). (8) 

The norm of the eigenvector is 

B^{X) = f-l = ^dj/{\-mif + l/cos2(7rA) 

i 

= d2s(A) + [l + dVr2(A)], (9) 

where (^) has been used. The first term comes from X^i^fj ^^^^ denoted as 
d^s(A), and the second term comes from X^n^n- unitary mixing matrix 
Uix in eqs. (|l|) and (||) is given by 

U,x^w,{\)lB{\). (10) 

For the ADDM model, whose mass matrix M' is not symmetrical, we need 
to compute the eigenvalues A^ and the eigenvectors i = (wi, • • • , w/, ui, W2, ■ • ■)"^ 
of AI'{M')'^ . The characteristic equation and the eigenvector components are 
once again given by (|^) and (^), if we set rrii — and n > 1. In addition, there 
are also / — 1 eigenvectors with A = 0, u„ = 0, and {wi, ■ ■ ■ ,Wf) orthogonal to 
idi,---:df). 
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4-1 Oscillation pattern 

We are now in a position to examine the consequence of this eigenstructure 
on the survival probabUhty Pair) given in (|). For the DDG model we will 
assume mi to be neither integers nor half integers, and all distinct. Since 
phases are absent in these models, only the real part contributes, hence 

ruiT) = l- J2iU^yU^x"fam^iAX^T/2). (11) 

A">A' 

4-2 Weak- coupling limit 

When di = 0, the matrix (||) is diagonal, with brane eigenvalues and bulk 
eigenvalues n S Z. For couplings d much smaller than the separation between 
any of these free eigenvalues, the eigenvalues are shifted very little so it is 
convenient to label them by the unshifted ones: ~ and A„ ~ n. The 
eigenvectors and their norms can be computed using perturbation theory from 
(§ to ® , to yield 

U^x, ^ 1 

UiXj ~ d^eiejlri{mj){mj - rm) 

U,x^^d^r{K)l{n-mi). (12) 

Assuming r(A„) and Tj{mj) to be of order unity, it folows that the small ele- 
ments UiXj and Uix^ are of order d"^ , and inversely proportional to the distance 
between the unshifted eigenvalues and nii. 

Appling this result to (pd]), we can approximate it by 

V.,ir)^l-J2ulsin'iAX^r/2). (13) 

A 

Thus the high 'frequency' components with large AA^ ~ A^ — mf are weakened 
by the factor 1/(A — mi)^. 

These results, derived for the DDG model, is also valid for the ADDM 
model with minor changes. Since / — 1 of the brane neutrinos decoupled from 
the rest of the ADDM model and can be taken into account easily, let us just 
discuss the case / = 1. In the absence of coupling, the eigenvalue is Ab = 
for the brane and A„ = n = 1, 2, • • • for the bulk. For weak coupling. A;, ~ d, 
Uixt ^ 1, and Uix„ ~ d^/n^. 
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4-3 Intermediate coupling 

In this case many KK modes are excited. The oscillation pattern gets very wig- 
gly and very complicated, as illustrated in Fig. 2 from a numerical calculation 
taken from Ref . [8] . It seems unlikely that this is the situation phenomenolog- 
ically unless we are so unlucky to have missed all these 'irregular' patterns so 
far. 

4-4 Strong coupling 

One might think that if intermediate coupling is difhcult to analyse, strong 
coupling will be impossible. This turns out not to be the case because the 
widths of the eigenstates become so large that they overlap with one another 
to eliminate most of the wiggly and complicated structure seen in intermedi- 
ate couplings. In fact, perhaps somewhat surprisingly, the problem becomes 
exactly soluble □ even for the DDG model. 

To explain how that comes about let us first examine the eigenvalues in the 
strong coupling limit {d ^ mj, 1). For any d > 0, we see from Fig. 1 that r(A) 
is a decreasing function of A, with poles at to^, but tan(7rA) is an increasing 
function of A, with poles at half integers. Consequently, there is one and only 
one solution of in any interval bounded by a neighboring pair of poles. As 
d is increased, the solution slides towards bigger magnitude of tan(7rA), i.e., A 
increases if r(A) > and decreases if r(A) < 0. When d oo, an eigenvalue 
must end up at the boundary of the interval, or at a zero of r(A). It is not 
difficult to show that the eigenvalues then consist of all half integers, plus the 
/ — 1 zeros of r{X). We denote the former by A^ = n + ^ (n G Z), and the 
latter by A^ {2 < £ < f). AJj will be referred to as the regular eigenvalues, and 
will be referred to as the isolated eigenvalues. Note that for large d, the 
brane and bulk neutrinos are quite thoroughly mixed up so we can no longer 
tell whether an eigenvector is more brane-like or more bulk-like. In particular, 
if we trace an isolated eigenvalue continuously back to d = 0, depending on 
the precise values of and to^ , sometimes it ends up as a brane eigenvalue nii , 
and sometimes it ends up as a bulk eigenvalue n. 

Next let us look at B^{\) and Uix. They behave differently for the two 
kinds of eigenvalues, because tan(7rA) is infinite for the regular eigenvalues and 
finite for the isolated eigenvalues. For this reason, it is more convenient to write 

for the regular eigenvalues in the form B^{\) — d'^s{\) + 1 + d"^ tan^(7rA), 
from which we see that B^{\) = 0{d^) and Uj\ is finite in the large-d limit. 
For regular eigenvalues, r(A) is non-zero, so we see from (||) that B{\) is of 
order d^. This means that Uj\ vanishes in the limit for every A of order unity. 
However, it does not mean that the sum in (||) over all regular eigenvalues is 
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zero, because there is an infinite number of regular eigenvalues. In fact, for 
large A of order of (P , we have r(A) ~ 1/A and s(A) ~ 1/A^, so B^{X) ~ 
1 + + 7r^d^)/A^ = 1 + K^/X"^ is finite for fixed d?/\ in the strong couphng 
limit. As a result, the sum in (H) over the regular eigenvalues can be replaced 
by the sum between n — ±N and ±00, for any finite number N. We will 
choose TV so that A'^ 3> 1,'Tii, then this sum can in turn be replaced by an 
integral, yielding 

00 



n— — oo 



00 ^2g-.A^r 



d\- 



A2 + ^2 

= e,e,5(ifV), (14) 

where 

g{x) ^ - / du^—- (15) 

begins at g(0) = 1. Its absolute value decreases monotonically in x, approach- 
ing g{x) ~ (1 — i)/^j2TTx for large x. 

This leaves the sum over isolated eigenvalues. Since Ui\ is finite for these 
eigenvalues, the transition in the strong coupling limit is amplitude 



s{\t){\i - mi){\e - rrij) 
V;,V,MK^r) + E V:,V,,e-^^'^^ (16) 



e=2 



It can be verified that the / x / matrix V is unitary (actually real orthogonal), 
so that Aij{0) = Sij as it should. If the function g{K^T) were replaced by 
e~*-^i'^ (for some Ai), then (^6|) would just be the usual transition amplitude 
of / fermions, with a mixing matrix Vij and no bulk neutrino. In that case 
unitarity of V guarantees the conservation of (brane) neutrino probability, 
namely, 7^j(r) = X]f=i^y(''') = 1 ^O'' With the presence of g{K'^T) in 

(H), there is a leakage into the bulk so the total brane neutrino probability 
decreases with 'time', as in 

7',(r) = l-|e,f{l-|.9(i^V)p}. (17) 
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This reaches an asymptotic amount 1 — jcj-p for large t. Since J^j kiP — 1; 
the average leakage per branc channel is 1//. 



5 Phenomenology 

According to the Super-Kamiokande analysis of their data0, the mixing angle 
for atmospheric neutrino is large, and the mixing with sterile neutrino is small. 
Their data on solar neutrino favor the large angle MSW solution LMA, without 
much mixing with a sterile neutrino. However, this conclusion on solar neutrino 
depends somewhat how the data are analysed. The small angle (SMA) and the 
vacuum oscillation (VAC) solutions, with or without a large amount of sterile 
neutrino present, are nat,completely ruled outM, even with the recent charge 
current data from SNOlla taken into account. 

There have-been several attempts to understand the data from bulk neu- 
trino models B'Ej. Unfortunately the data pose considerable difficulties on the 
simple DDG/ADDM model, which contains only one bulk neutrino and has 
no direct mixing between the active neutrinos. In the weak coupling limit 
when the mixing is small, there is no way this simple model can be applied 
to the atmospheric neutrino problem. Jvlotivated by the 2+2 solution intro- 
duced to explain also the LSND datatJ, in which {v^^v^) and {Vf.,Vs) form 
nearly degenerate doublets separated by about 1 eV, and the solar neutrino 
oscillates mostly to the sterile neutrino, attempts have been made to use the 
DDG/AGGM model to explain the solar neutrino problem. However, the 2+2 
solution contains one sterile neutrino Vg, but the DDG/ADDM model contains 
an infinite tower of them. With more sterile neutrinos present the oscillation 
pattern becomes more complicated, and this can be utilized to cure the wrong 
energy dependence in the Super-K data caused by a single Ug. Two groups 
have succeeded in fitting the solar neutrino data, befora.the SNO data be- 
came available anyway, either with a / = 1 DDG modelLJ,|jQr by introducing 
a mass for the bulk neutrino in the / = 1 ADDM model 113. To be able to 
explain also the atmospheric and perhaps the LSND data, one must mod- 
ify the DDG/ADDM model, by introducing direct mixing between the active 
neutrinos, by introducing three bulk neutrinos, or by some other means. 

One might think that the large mixing of atmospheric neutrino, and per- 
haps also of the solar neutrino, could possibly be explained by the DDG/ADDM 
model in the presence of a strong coupling. However, strong coupling also di- 
verts a considerable amount of active neutrino flux into sterile neutrinos (the 
first term within the square brackets in (16)), which is forbidden by the data. 
Nevertheless, if the diversion comes mostly from an incoming Vr beam, for 
which there is no data, it is conceivable that it might still work. In other 



9 



words, we might attempt to make e^- in (|T^) close to 1 and ee,e^ small, in 
which case there is not much leakage of the solar and atmospheric neutrino 
fluxes into the sterile channels. Unfortunately this does not work for the / = 3 
DDG model, because with the vanishing of the first term of ( [l6| ) for r ^ ^/K^, 
there remains only / — 1 = 2 terms in ( p^ and only one square-mass differ- 
ence, impossible to explain both the atomospheric and the solar oscillations. 
Whether there is a way to make this work for / > 4 has not yet been analysed. 
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Figure 1: Graphical solutions of the characteristic equation (|^). The left hand side of (|6|) 
is represented by the solid curve; the right hand side is represented by the dashed curves, 
with (mi, 7712, ms) = (—1.8,-0.2,0.3). The one with the short dash has the larger coupling 
constant d than the one with the long dash. Note that the solution A moves up and right 
along the solid curve for increasing d when the value on both sides of ^ is positive, and 
moves down and left along the solid curve when the value is negative. 
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Figure 2: Taken from Fig. 5 of Ref. [8] to show the wiggly pattern in the L/E dependence, 
with a moderate couphng to the infinite tower of KK neutrinos. The boxes are the Super- 
Kamiokande and K2K results. 
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